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Abstract 

This work discusses the homogenization analysis for diffusion processes on scale-free metric graphs, using 
weak variational formulations. The oscillations of the diffusion coefficient along the edges of a metric graph 
induce internal singularities in the global system which, together with the high complexity of large networks 
constitute significant difficulties in the direct analysis of the problem. At the same time, these facts also 
suggest homogenization as a viable approach for modeling the global behavior of the problem. To that 
end, we study the asymptotic behavior of a sequence of boundary problems defined on a nested collection 
of metric graphs. This paper presents the weak variational formulation of the problems, the convergence 
analysis of the solutions and some numerical experiments. 
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1. Introduction 

A scale-free network is a large graph with power law degree distribution, i.e., IP[deg(n) = fc] ~ k~^ 
where 7 is a fixed positive constant. Equivalently, the probability of finding a vertex of degree k, decays 
as a power-law of the degree value. Power-law distributed networks are of noticeable interest because they 
have been frequently observed in very different fields such as the World Wide Web, business networks, 
neuroscience, genetics, economics, etc. The current research on scale-free networks is mainly focused in 
three aspects: first, generation models (see mm, second, solid evidence detection of networks with power- 
law degree distribution (see 0IH El E] ) ■ The third and final aspect studies the extent to which the power-law 
distribution relates with other structural properties of the network, such as self-organization (see mi); this 
is subject of intense debate, see [7] for a comprehensive survey on the matter. 

The present work studies scale-free networks from a very different perspective. Its main goal is to 
introduce a homogenization process on the network, aimed reduce the original order of complexity but 
preserving the essential features (see Figure^. In this way, the “homogenized” or “upscaled” network is 
reliable for data analysis while, at the same time, involves lower computational costs and lower numerical 
instability. Additionally, the homogenization process derives a neater and more structural picture of the 
starting network since unnecessary complexity is replaced by the average asymptotic behavior of large 
data. The phenomenon known as “The Aggregation Problem” in economics is an example of how this 
type of reasoning is implicitly applied in modeling the global behavior of large networks (see [9]). Usually, 
homogenization techniques require some assumptions of periodicity of the singularities or the coefficients of 


*This material is based upon work supported by project HERMES 27798 from Universidad Nacional de Colombia, Sede 
Medellin. 

* Corresponding Author 

Email address: fainoralesj@unal.edu.co (Fernando A Morales) 


Preprint submitted to arXiv 


May 31, 2016 





the system (see nnmn), in turn this case demands averaging hypotheses in the Cesaro sense. The resulting 
network has the desired features because of two characteristic properties of scale-free networks. On one 
hand, they resemble star-like graphs (see |12jl. on the other hand, they have a “communication kernel” 
carrying most of the network traffic (see [T5] l. 


This paper, for the sake of clarity, restricts the analysis to the asymptotic behavior of diffusion processes 
on stared metric graphs (see Definition^emd below). However, while most of the models in the 

preexistent literature are concerned with the strong forms of differential equations (see m for a general 
survey and for the stochastic modeling of advection-diffusion on networks), here we use the variational 
formulation approach, which is a very useful tool for upscaling analysis. More specifically, we introduce the 
the pseudo-discrete analogous of the classical stationary diffusion problem 


-V-{K\/p)=f inH, 
p = 0 on dV ,, 


( 1 ) 


where K is the diffusion coefficient (see Definition ^a.nd Equation (10) below). Due to the variational for¬ 
mulation it will be possible to attain a-priori estimates for a sequence of solutions, an asymptotic variational 
form of the problem and the computation of effective coefficients. Finally, from the technique, it will be 
clear how to apply the method to scale-free metric networks in general. 


Throughout the exposition the terms “homogenized”, “upscaled” and “averaged” have the same meaning 
and we use them indistinctly. The paper is organized as follows, in Section^tde necessary background is 
introduced for L^, H^-type spaces on metric graphs as well as the strong form and the weak variational 
form, together with its well-posedness analysis. Also a quick review of equidistributed sequences and Weyl’s 
Theorem is included to be used mostly in the numerical examples. In Section we introduce a geometric 
setting and a sequence of problems for its asymptotic analysis, a-priori estimates are presented and a type 
of convergence for the solutions. In Section under mild hypotheses of Cesaro convergence for the forcing 
terms, the existence and characterization of a limiting or homogenized problem are shown. Finally, Section 
[^is reserved for the numerical examples and Section holds the conclusions. 




(b) Homogenized Network. 


Figure 1: Figure (a) depicts a scale-free network. Figure (b) depicts a homogenization of the original network. 
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2. Preliminaries 


2.1. Metric Graphs and Function Spaces 

We begin this section recalling some facts for embeddings of graphs. 

Definition 1. A graph G = (V^E) is said to be embeddable in if it can be drawn in so that its 
edges intersect only at their ends. A graph is said to be planar if it can be embedded in the plane. 

It is a well-known fact that any simple graph can be embedded in IR,^ or IR^ (depending whether it is planar 
or not) in a way that its edges are drawn with straight lines; see |16j for planar graphs and |17j for non-planar 
graphs. In the following it will always be assumed that the graph is already embedded in IR^ or IR^. 

Definition 2. Let G = {V, E) be a graph embedded in IR^ or IR^, depending on the case. 

(i) The graph G is said to be a metric graph if each edge e € E is assigned a positive length G (0, oo). 

(ii) The graph G is said to be locally finite if deg(u) < -l-oo for all v € V. 

(iii) If the graph G is metric, the boundary of the graph is defined by the set of vertices of degree one. 
The set will also be denominated as the set of boundary vertices and denoted by 

dV {v & V : deg(u) = 1}. (2) 

(iv) Given a metric graph we define its natural domain by 

U int(e). (3) 

ee E 


Definition 3. Let G = (V, E) be metric graph we define the following associated Hilbert spaces 
(i) The space of square integrable functions, or mass space is defined by 


LHG) 0 L"(e), 

ee E 

endowed with its natural inner product 

{f:9)L^G) '^ [ f 9- 


(4a) 


(4b) 


(ii) The energy space of functions is defined by 

iji(G) {/ e 0 H\e) : Jim /(a:) = Jim /(a;), 

E X Ge X £ (T 

for all vertices u G V and for all edges e,a incident on u|. (5a) 

In the sequel f{v) lim{/(a:) : x ^ v, x G e}, with e G E any edge incident on v. We endow the 
space with its natural inner product 

{f,9)m(G) 51 [ ^9 + (5b) 

Here de denotes the derivative along the edge e G E. 
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(iii) The space Hq{G) is defined by 


ffo(G) =' { / G H\G) : f{v) = 0, for all v&dV], 
endowed with the standard inner product ( |5b[ ). 

Remark 1. Let G be a metric graph 


( 6 ) 


(i) Notice that the definition of de is ambiguous in the expression (5b), such ambiguity will cause no 
problems since the bilinear structure of the inner product is indifferent to the choice of direction 


(q,r) deqd^r = ( - d^q) {-dev). 


(ii) Whenever there is need to specify the direction of the derivate, we write de,v to indicate the direction 
pointing from the interior of the edge e towards the vertex v on one of its extremes. 

(iii) Notice that if the metric graph G is connected, then the Poincare inequality holds and the inner product 


(iv) 




( 7 ) 


eG 


is equivalent to the standard one (5bI in the space Hl{G). 


Observe that the condition of agreement of a function / G (G) on the vertices of the graph G does not 
necessarily imply continuity as a function / : Oq —>■ IR. For if the degree of a vertex u G F is infinite and 
the function is continuous on v, then it follows that the convergence f(v) lim{/(a;) : a; —>■ u, x G e} 
is uniform for all the edges e incident on v. Such a condition can not be derived from the norm induced 
by the inner product (5b), although the function /lint(e) is continuous for all e G E. 


Definition 4. Let G„ = (F„,i?„) be a sequence of graphs. 

(i) The sequence {G„ : n G N} is said to be increasing if Vn C Vn+i and Q E^+i for all n G N. 

(ii) Given an increasing sequence of graphs {G„ : n G N}, we define the limit graph G = {V,E) in the 
natural way i.e.. 


U f. 

n G ]N 


y E, 

n G ]N 


In analogy with monotone sequences of sets we adopt the notation 

Gli' U G„. 

n G IN 


2.2. The Strong and Weak Forms of the Stationary Diffusion Problem on Graphs 

Definition 5. Let G = (F, E) be a locally finite metric graph, F G L^{G) and h : V — dV — >■ IR, define the 
following diffusion problem 

X -de{K9ep)le = X in Hq. (9a) 

e G S ee E 

Where K G L°°{nG) is a nonnegative diffusion coefficient. We endow the problem with normal stress 
continuity conditions 

Jim^p(x) = p{v) for all p G F — dV, (9b) 

X G e 


4 





and normal flux balance conditions 

^ =0 for all u G F — dV. (9c) 

e^E a; G e 
e incident on v 

Here de,v denotes the derivative along the edge e pointing away from the vertex v. Finally, we declare 
homogeneous Dirichlet boundary conditions 

p(v) = 0 for all V G dV. (9d) 

A weak variational formulation of this problem is given by 

peH^iG): E Hv)q{v) yqeH^iG). (10) 

eG-E'^® eeE''^ veV-dV 


For the sake of completeness we present the following standard result. 

Proposition 1. Let G = {V,E) be a locally finite connected metric graph such that dV ^ 0 and let K G 
L°°{LIg) be a diffusion coefficient such that K{x) > ck > 0 almost everywhere in LIq. Then the problem 
(101 is well-posed. 


Proof. Clearly the functionals on the right hand side of Problem (10) are linear and continuous, as well as 


the bilinear form b(p, q) g e /e ^^eP deq of the left hand side. Additionally, 


^ J >ckY1 


e G E ‘ 


e G E ‘ 


\deP\^ >cJ2 Mn^e) = c|b||ffi(G)- 

e G E 


The first inequality above holds due to the conditions on K. The second inequality hods due to the Dirichlet 
homogeneous boundary conditions and the connectedness of the graph G, which permits the Poincare 
inequality on the space Hq{G) as discussed in Remark (|^ above. Therefore, due to the Lax-Milgram 
Theorem, the Problem (10) is well-posed. □ 


2.3. Equidistributed Sequences and Weyl’s Theorem 

The brief review of equidistributed sequences and Weyl’s theorem of this section will be applied, almost 
exclusively in the numerical examples below, see Section For a complete exposition on equidistributed 
sequences and Weyl’s Theorem see ng. 

Definition 6 . A sequence {0„ : n G N} is called equidistributed on an interval [a, b] if for each subinterval 
[c, d] C [a,b] it holds that: 

#{i : Oj G [c,d], l<i<n} _ d-c 
n->oo n b — a 

Theorem 2 (Weyl’s Theorem). Let {d„ : n G N} be a sequence on [a,b], the following conditions are 
equivalent: 

(i) The sequence : n G N} is equidistributed in [a, &]. 


(ii) For every Riemann integrable function / : [a, &] —>■ C 


lim - ^ f{e,) = 
n—>-oo n ^^ 0 — a 

2=1 


f{0) de. 


Definition 7. Let H = i?(0, 1) C and let / : H > IR be such that its restriction to every sphere dB(0, p) 
with 0 < p < 1 is Riemann integrable. Then, we define its angular average by the average value of / 
along the sphere 9(R(0,p)), i.e.. 


m,[/] : [0,1)^]R, 


(0 7 ^ / f{tcos0, tsind)) dO. 

Stt ./n 


( 12 ) 
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3. The Sequence of Problems 

In this section we analyze the behavior of the solutions {p" : n G N} of a family of well-posed problems 

on an very particular increasing sequence of graphs {Gn : n G N}, depicted in 

3.1. Geometric Setting and the n-Stage Problem 

In the following we denote by 17, 5'^ the unit disk and the unit sphere in respectively. The function 

F : VIq —7^ IR is such that T'la„ G LS‘{Vln) for all n G IN, {ft," : n G N} is a sequence of real numbers and 

the diffusion coefficient K G L°°{VIq) is such that K{-) > ck > 0 almost everywhere in 17^. 

Definition 8. Let {u„ : n > 1} be an equidistributed sequence in and vq 0 G IR^. 

(i) For each n G N define the graph G„ = (14, En) in the following way: 

v„ {Vn -.0 <i < n}, En {vqVi : 1 < z < n}. (13) 

(ii) For the increasing sequence of graphs {G„ : n G N} define the limit graph G Une jn described 

in Definition 

(iii) In the following we denote the natural domains corresponding to G, Gn by Qq ^^nd f7„ respectively. 

(iv) For any edge e G E we denote by Ve its boundary vertex and 9e G [0, 27r] the direction of the edge. 

(v) From now on, for each edge e = voVe and / : e ^ IR a function, it will be understood that its 
one-dimensional parametrization, is oriented from the central vertex vq to the boundary vertex Vg. 
Consequently the derivative dg equals dg^v^- 

(vi) For any given function / : flc —> IR (or f ■. ^ fR-) we denote by fg : (0, 1) —>■ IR, the real variable 

function fg(t) (/le)(tcosfte,f sinde) on the edges e £ E {oi e £ En respectively). 




Figure 2: Figure (a) depicts the stage 5 of the graph G. Figure (b) depicts a more general stage n of the graph G. 

Remark 2. From the following analysis, it will be clear that it is not necessary to assume that the sequence 
of vertices {u„ : n G N} of the graph is equidistributed or that the vertices are in or even that the graph is 
embedded in IR^. We adopt these assumptions, mainly to facilitate a geometric visualization of the setting. 
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From now on it will be assumed that {G„ : n € N} is the increasing sequence of graphs, with G its limit 
graph, as in the Definition^ above. Next, we define the family of problems to be studied, for each n S N 
consider the well-posed problem 

p^GH^iGn)-. Y. fKdeP^deq=Y ! F q + q{vo) , \JqGHl{Gn). (14) 

We are to analyze the asymptotic behavior of the sequence of solutions {p” : n € N}. One of the main 
challenges is that the elements of the sequence are not defined on the same global space. The fact that p"(0) 
may not be zero makes impossible to extend this function to Hq(G) directly, however it will play a central 
role in the asymptotic analysis of the problem. 


3.2. Estimates and Edgewise Convergence Statements 

In this section we get estimates for the sequence of solutions, several steps have to be made as it is not 
direct to attain them. We start introducing conditions to be assumed from now on 


def 


Hypothesis 1. (i) The forcing term F is defined in the whole domain, i.e. F : ^Iq —> IR and M — 

SUPeeB ||-f^||L2(e) < +00. 

(ii) The sequence | — /i" : n € N } is bounded. 

^ n ^ 

(hi) The permeability coefficient satisfies that K G L°°(Qg), infxeQo K(x) = ck > 0 and iFle = K{e) 
i.e., it is constant along each edge e G E. 

Remark 3. Notice that the Hypothesis [^(ii) states that the balance of normal flux on the central vertex 
is of order 0(n) i.e., it scales with the number of incident edges. 

Lemma 3. Under the Hypothesis^ the following facts hold 

(i) The sequence {p"(0) : n G N} C ]R is bounded. 

(ii) Let e G E be an edge of the graph G then, the sequence {deP^(0) : e G En} C ]R zs bounded. Moreover, 
there exists Mq such that |9eP”(0)| < Mq for all e G E and n G N such that e G E^. 

(Hi) Suppose that the sequences { — lo i^~^)Fe{t)dt :nG]N},{ — /i”:nGN} and {— J2eeE„ Fi{e) : 


n G N} are convergent, then the following limits are satisfied 


lim p"(0) = lim 

n—¥oo n—>'OC 

For any fixed edge e G E holds 


F SeGB„ Jo - ^)Fe(t)dt - 

n SeGB„ Fi{e) 


pl 

lim K{e)deP^{0) = L{e) / (t — l)Fe(t)dt — K{e) lim 

n—^oo Jq n—^oc 


i EaeE„ lo i* - ^)Fa{t)dt - 


nT.a&E„K{a) 

Moreover, the convergence is uniform in the following sense 

Ve>0 3N gTN such that, ifn > N and e G En, then \K{e) deP^iO) — L{e)\ < e. 


(15a) 


(15b) 


(15c) 


Proof. (i) Let q G Rq(G„) be the function such that q{0) = —1 and qe{t) = t — 1 for all e G i?„. Test 


(141 with q, this yields 


g(0) = 


e G 


e G 


Fq + K^q{b). 


(16) 
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Computing and doing some estimates we get 


#Er.CK\p^m<^ E \\F\\L^ie) + \h 
eeE„ 


Hence 


This proves the first part. 


b"(0)| < —M+ — 
Ck Ck 


(ii) Let e G E he a, fixed edge, let n G N be such that e G En and let p" be the solution to Problem (14). 
Let q G iLg(G'„) be as in the previous part and test (14) to get 

-iL(e)p"(0)g(0)+ E d,p^d,q = K{e)p'-{Q)- E 9 + <z(0) E 

^ ^ ■H' J (7 _^T71 


(T G E. 
o^e 


(7 G En 
(7^e 


<7 G En 
(7^e 


= j FqP Y. [ Fq + h'^q{0). 

J e ^ ^ rr J (7 


(7 GE, 
(7^e 


In the expression above, integration by parts was applied to each summand a ^ e oi the left hand side, 
in order to get the second equality. Now, recalling that J2eeE„ F dep'^{h) = and that —Klgd'^p'^ = 
Ele for each e G En the equality above reduces to 


K(e)p"(0)= [ it-l)F,{t)dt-Kie)deP’^{0). 
Jo 


Hence, 


\deP^m < b”(0)| + — \\F\\me) < —M + — 

CK CK CK 


(17) 

(18) 


Choosing Mg > 0 large enough, the result follows. 


(hi) Let q G Hl{Gn) be as in the previous part, testing (14) with it yields the equality (16) which is 
equivalent to 


P 


e G 


ee E„ 


Now, letting n —>■ oo the equality (15a) follows because the hypothesis K{e) > ck for all e G i? implies 
that — g K{e) >CK>h. For the convergence of {£?eP"(0) : e G ifn}, let n —)■ oo in the expression 


n J2^eE^ fo i)Ea(t)dt 


to get the equality ( |15b[ ). For the uniform convergence observe that the identity ( |17[ ) yields 

- ^(e)p”(0) 


|iG(e)9ep"(0) — L(e)| = K{e) lim 


<l!^lh 


lim 

n—^oo 


lEaeE^F{a) 

^ EaeE^ fo(* - i)F^(t)dt - 

^EaeE^F(a) 


-pEo) 


Finally, choose iV G N such that the right hand side of the expression above is less than e > 0 for all 
n > N then, the term of the left hand side is also dominated by e > 0 for all n > and e G En- □ 

Remark 4. It is clear that in Lemma^paii (hi), it suffices to require the mere existence of the limit 

E<tg£;„ fo(* - l)F<T(t)dt - 


lim 

n—>-oo 


EaGEn F{<^) 


in order to attain the same conclusion. However, the hypotheses of (iii) are necessary to identify the 
asymptotic problem and compute the effective coefficients. 



















Theorem 4. Let F, {/i" : n € N} and K verify the Hypothesis^ as in Lemma^ then 
(i) There exists a constant Mi such that < Mi for all e G E and n G IN such that e G E^- 

(ii) For each e G E there exists G H^{e) such that ||p"le ~ ^ 0- Moreover, this 

convergence is uniform in the following sense 

Ve > 0 3N G N such that, ifn > N and e G E„, then \\deP^ — deP^^'^Wn^ie) < e- (19) 


(Hi) The function p : Lie IR- given by pig is well-defined and it will be referred to as the limit 

function. 


Proof. (i) Fix e G E and let n € IN be such that e G E^. Since p” is the solution of Problem (14) it 
follows that —K{e) dlp^ = Fig € F^(e) for all e S in particular p”lg S F^(e) with ||£?eP"||L 2 (g) < 

— ll^lU^(e) < — M. On the other hand, since c?eP" le is absolutely continuous, the fundamental 
Ck Ck 

theorem of calculus applies, hence 5gp”(x) = 9gp"(0) + d'^p'f{t) dt = 0gp”(O) + Ff,{t) dt for all 

X G e. Therefore, 

|agp"(a;)p = 2|0gp"(O)p + 2a;||F||i.(g) < 2M^ + 2M^. 

Where Mg is the global boun d found in Lemma [^(ii) above. Integrating along the edge e gives 
\\9eP’"\\L2(e) < \/2(M^ + M2). Next, given that p^{v) = 0 for all v G E^, repeating the previous 

argument yields |lp"||L 2 (g) < sJ2{Mq + M^). Finally, since ||9gp”||L2(e) < — M, the result follows 

^K. 

for any Mi satisfying 

Ml > AMl + (4 + -\r)M^. 

V c^/ 

(ii) Fix e G E, due to the previous part the sequence {p"le : e G F„} is bounded in H^{e), then there 
exists G H‘^{e) and a subsequence {uk : fc G N} such that 

p^k - y p{e) -^^eakly iniL^(e) and strongly iniL^(e). 


Let (fi G H^{e) such that equals zero on the boundary vertex of e. Let q be the function in Hl(Gn) 
such that pe = <p and qcr{t) = ‘p(0)(l — t) is linear for all a G En — {e}. Test Problem (14) with this 
function to get 

[ K{e)d,p^d,qF Y. ^ f f FQ+ Y [ Fq + h^q{0). 


<y ^En 
a^e 


(7 G Ert 
o^e 


Integrating by parts the second summand of the left hand side yields 


iF(e)agp" Sgp - Y^ dlp^q-g,{0) Y K d^p^O) = Fip+ Y / Fg, + h^q{0). 


<7 ^ En 
(7^e 


<7 G En 

(7^e 


(7 G En 
(7^e 


Since p" is a solution of the problem, the above reduces to 

J K{e)dep'^ dep + K{e)dep'^{0) (p{0) = j Fq. 


( 20 ) 


Equality (20) holds for all n G N, in particular it holds for the convergent subsequence {uk : fc G N}, 
taking limit on this sequence and recalling (|15b ), we have 


J^K{e)deP^''^ deP = J Fp - L{e)p{0). 


(21) 












The statement holds for all ip G H^{e) vanishing at Ve, the boundary vertex of e. Define the space 
H{e) {p G H^{e) : p{ve) = 0} and consider the problem 


uGH{e): J K{e)deudeP = J Fp — L{e) p{0), \/pGH{e). 


( 22 ) 


Due to the Lax-Milgram Theorem the problem above is well-posed, additionally it is clear that G 
H{e), therefore it is the unique solution to the Problem (22) above. Now, recall that {p"le : e G F^} 
is bounded in F[^{e) and that the previous reasoning applies for every strongly 77^(e)-convergent 
subsequence, therefore its limit is the unique solution to Problem (22). Consequently, due to Rellich- 
Kondrachov, it follows that the whole sequence converges strongly. Next, for the uniform convergence 
test both statements (20) and (21) with (p^le —and subtract them to get 

Cifll^eP” - OeP^®)||ffi(e) < K {e) f | ^ep" - = (L(e) - R'(e)aep" (0)) (p” (0) - p^^^ (0)) 

J e 

< \L{e) - K{e)deP^(0)\ H^eP” - aeP^^^lkqe)- 


The above yields 


llOeP” - < — \L{e) - (c) ^ep” (0) I. 


CK 


Now, the uniform convergence (19) follows from the Statement (15c), which concludes the second part, 
(hi) Since p*^®^(0) = lim p”(0) for all e S i? then, the limit function p is well-defined and the proof is 

n—^oo 

complete. □ 


4. The Homogenized Problem: a Cesaro Average Approach 


In this section we study the asymptotic properties of the global behavior of the solutions {p" : n € N}. 
It will be seen that such analysis must be done for certain type of “Cesaro averages” of the solutions. This 
is observed by the techniques and the hypotheses of Lemma which are necessary to conclude the local 
convergence of {p"le '■ e G En}. Additionally, the type of estimates and the numerical experiments suggest 
this physical magnitude as the most significant for global behavior analysis and upscaling purposes. We 
start introducing some necessary hypotheses. 

Hypothesis 2. Suppose that F, {/i" : n G IN} and K verify Hypothesis ^and, additionally 

(i) The diffusion coefficient K : Lie (OjOo) has finite range. Moreover, if K{E) = {Ki : 1 < z < /} and 
B, {eG E: K{e) = AJ, then 


e G E„nBi 


#(A„ n Bi) 

I\q 

n 


-G Si Ki. 


(23) 


With Si > 0 for all 1 < z < J and such that ^ Sj = 1. 


(ii) The forcing term E satisfies that 

1 


#(A„ n Bi) 


E 


Fi , 


V 1 < z < /. 


(24) 


eG E„nBi 


Where Fi G T^(0,1) and the sense of convergence is pointwise almost everywhere. 

(hi) The sequence | — /z" : rz G N| is convergent with h = lim — h". 

rz ^ n-zoo ri 
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Remark 5. (i) Notice that if (i) and (ii) in Hypothesis ^are satisfied, then 




#{En n Bi) 


1 


n #iEn n Bi) 


E 


s^Ei. 


eGB„ i=l \ eeE^nBi 

Hence, the sequence {Ff, : e G E} is Cesaro convergent. 


i=l 


(ii) A familiar context for the required convergence statement (24) in Hypothesis^above is the following. 
Let A be a continuous and bounded function defined on the whole disk H and suppose that for each 
1 < i < /, the s^uence of vertices : n S N, VuVq G Bi} is equidistributed on S^. Then, due to 

for any fixed t G (0,1) it holds that ^ p ^ J2ee E„nBi ^e{t) -S' ing)/] i.e, 


Weyl’s Theorem 
the angular average introduced in Definition 


fi{E„ n Bi) 


4-1. Estimates and Cesaro Convergence Statements 

Lemma 5. Let F, {/i" : n G N} and K verify Hypothesis^ then 

(i) The sequence { — ^ p'f : n G is bounded in H^(0,1). 


(ii) The sequence 


■ #{En n Bi 


Y p” : n e NI 


(25) 


e G E„nBi 


is bounded in H^{0, 1) for all i G {!,..., /}. 

Proof. (i) Test Problem ( |l4| ) with p” to get 

CK liaeP"lli2(,) < 5] fK\dp-\^= Y [Fep- + h-p-M 


e G 


eG 


e G 


< 


( E ll^’lli^(e))'^'( E ll^^”lli^(e))''E|/l’^||p”(uo)|. 

e G e G 

Since p"(ve) = 0 for all e G £!„, then |b"|U 2 (g) < \\dp’^\\L 2 (e) and ||p”||ffi(e) < ||0p"|U2(g). Hence, 

dividing the above expression over n gives 

E E ii^«i-(.))'''(s E iij-”iit.(.))‘'A: 


1/2 


1 

n 


eG 


< 


e G B- 

M /I 


e G B„ 




^£(i E 


CK '^n 


e G B„ 


Here C > 0 is a generic constant independent from n G IN. In the second line of the expression 
above, we used that M = supgg^^ , \\F\\L^(e) < + 00 , /i” : n G N} are bounded and that {p"(uo) : 

n G N} is convergent (therefore bounded) as stated in Lemma [^(i). Hence, the sequence Xn 
(^EeGB„ is such that < 2 ^ Xn + C for all n G IN, where the constants are all 

non-negative. Then {a;„ : n G IN} must be bounded, but this implies 


^ ^ - E IbeliBHO,!) = " E lb”II^Te) < (“ E 

n Bdo.i) n ^ ' n ^ ’ \n 


eGB„ ' ^ eGB„ 

Finally, recalling the estimate 




ee E„ 


e G B,i 


1/2 


CK 

n 

e G B„ 

the result follows. 




M, 


e G B„ 


e G B„ 
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(ii) Fix i e {1, 2,..., then 


1 


#{En n Bi) 


E 

eG E„nBi 


e t -C^Ti 


On the right hand side term of the expression above, the first factor is bounded due toHypothesis 
(iii), while the boundedness of the second factor was shown in the previous part. Therefore, the result 
follows. □ 

Before presenting the limit problem we introduce some necessary definitions and notation 
Definition 9. Let K verify Hypothesis^and I = =f^K{E) then 

(i) For all 1 < J < / define wt (cos(^ i), sin(^ f)) G S^, wq = 0 and V {wi : 0 < i < I}. 


(ii) For all 1 < J < / define the edges Ci woWi and £ {at ■ 1 < i < !}■ 


def 


(iii) Define the upscaled graph by Q 

(iv) For any (p G Hq{Q) and n G N denote by Tnp G iLQ(G„), the function such that (Tnp)le agrees with 
pl-ai whenever e & Bi. This is summarized in the expression 




e^E„ i = leG E„r\Bi 

In the sequel we refer to Tnp as the ffg (G„)-embedding of p. 


(v) In the following, for any 1 < I < J, we adopt the notation di da^ ■ Similarly, for any given 

function / : Og —>• IR and edge ai G £ we denote by fi : (0,1) IR, the real variable function 

Mt) (,/l<T.)(lcos(^I),Isin(^I)). 

Theorem 6. Let F, {/i" : n G N} and K verify Hypothesis^ then 
(i) The following problem is well-posed. 

I . I 


pG Hq{Q) : 


/ s^Ki d^p diq^^ / s^ Fi Pi + h(0) q(0) 


VqGHo^(g). 


(26) 


In the sequel, we refer to Problem (261 as the upscaled or the homogenized problem and its solution 
p as the upscaled or the homogenized solution indistinctly. 

(ii) The sequence of solutions {p" : n G N} satisfy 

1 


#(En n Bi) 


E Pe-Pi 


e G 

(iii) The limit function p : LIq —^ IR satisfies 

1 

n R. i 

e G E„nBi 

1 


n-^oo 


GO, 


#(En n Bi) 


Pe - Pi 


n^oa 


> 0 , 


E “ E 

e £ En i — 1 


V 1 < i < /. 


(27) 


V 1 < I < /. 


^1(0,1) n-^-oo 


GO. 
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Proof. (i) It follows immediately from Proposiiion 
(ii) Let (p S Hq{Q) and let Tn(p be its iLg(G'„)-embedding. Notice the equalities 


1 

n 


/ KdeP: 9eTnP = Y. ^ Y / 


e G ® 


eG S„nBi 


_ #(£'n n Bj) 

Z_^ n 


IK 


1 


0 \#{En(^Bi) 


Y p") ’ 

e G E^nBi 


and 


-1: / f T,*. = j;: 


eG En 


2=1 


#iEn n B,) /■! / 1 


E 

e G E^nBi 


Now, testing Problem (141 with 


— TnV^, due to the previous observations gives 
n 


(28a) 


(28b) 


E 

2=1 


#{En n Bi) 


'^‘Yi 


0 ^#{EnPBi) 


Y p") 


e G Er^nBi 


#(-E'ra n Bj) 


2 = 1 


K 


0 V#(.B„nBj) 


E f)^ 


eG E„nBi 


-p{0). (29) 


Due to Lemma [^(ii) there must exist a subsequence {n^ : k G IN} and a collection : 1 < i < d} C 
iL^(0,1) such that 


#{Enk n Bi) 


Y 


-G ^i weakly in iL^(0,1) and strongly in iL^(0,1), for all 1 < * < /. 


eG E^.nBi 


On the other hand, due to Hypothesis^{u) the integrand of the right hand side in the identity (291 is 
convergent for all i G {1,..., Due to Hypothesis 2l(i) the sequences : n G N} are also 


convergent for all i G {1,...Then, taking the equality (29l for and letting k ^ oo gives 


Y^iK, d^^dpi = Y^^ Fipi + hp{0). 


2= 1 


2=1 


Notice that since p"'' G Hl{Gnk) then 
1 


#{En^ n Bi) 


Y p:M 0)= ..p E ^'eHO), V*,JG{1,.(30) 


eG E„,nBi 


In particular ^i(O) = ^j(O); consequently the function r] G Hq{Q) such that iji = ^i is well-defined. 
Moreover, the identity (30) above is equivalent to 


Y / 9iP dip = E / P + ■ 


2 = 1 ' 


2 = 1 ' 


Since the variational statement above holds for any p G Hl[Q) and t] G Hq{Q), due to the previous 
part it follows that r] = p. Finally, the whole sequence {p" : n G Nj satisfies ( [2^ because, for 
every subsequence : j G N} there exists yet another subsequence {p”^« : £ G N} satisfying the 
convergence Statement (|^. This concludes the second part. 


13 



















(iii) Both conclusions follow immediately from the previous part and the uniform convergence Statement 
(191 shown in Theorem □ 

Remark 6 (Probabilistic Flexibilities of the Results). Consider the following random variables 

(i) Let X : E' —)■ (0, oo) be a random variable of finite range {Ki : 1 < i < /} and such that IE[X = Kj\ = Si 
for all 1 < i < /. Notice that due to the Law of Large Numbers, with probability one it holds that 


- E 


Si Ki. 


eG E„nBi 


(ii) Let Y : E ^ 1) be a random variable such that supgg^ ||y(e)||i 2 (g) < +oo and such that 


#(s„ n Bi) 


E 




V 1 < i < /. 


(31) 


(32) 


e G E^nBi 


Therefore, the results of T/ieorem hold, when replacing E by X or E by F or when making both substi¬ 
tutions at the same time. 


5. The Examples 

In this section we present two types of numerical experiments. The first type are verification examples, 
supporting our homogenization conclusions for a problem whose asymptotic behavior is known exactly. The 
second type are of exploratory nature, in order to gain further understanding of the phenomenon’s upscaled 
behavior. The experiments are executed in a MATLAB code using the Finite Element Method (FEM); it is 
an adaptation of the code femld.m [TO] . 


5.1. General Setting 

For the sake of simplicity the vertices of the graph are given by vi (cos £, sin £) G S^, as it is known 

that {ve : £ S N} is equidistributed in S^ (see [IH]). The diffusion coefficient hits only two possible values 
one and two. Two types of coefficients will be analyzed, Kj_,Kp a deterministic and a probabilistic one 
respectively. They satisfy 


Kd-. no ^{1,2}, 

Kdivgvo) '=*^ 12 ’ 

£ = 0 mod 3, 

£ ^ 0 mod 3. 

(33a) 

Kp : no ^ {1,2} , 

nKp = 1] = ^, 

m[Kp = 2] = ^. 

(33b) 


In our experiments the asymptotic analysis is performed for Kp being a fixed realization of a random 
sequence of length 1000, generated with the binomial distribution 1/3, 2/3. Since #Kd{E) = ^Kp{E) = 2 
it follows that the upscaled graph Q has only three vertices and two edges namely Wi = (1,0), ^2 = (—1, 0), 
Wo = (0,0) and ui = wiWq, <J 2 = W 2 W 0 . Also, define the domains 


= U e N , K{v(,Vo) = 1 } , = U ^ G N , K{veVo) = 2} . 

Where K = Kd or K = Kp depending on the probabilistic or deterministic context. Additionally, we define 


def 

Pi = 


1 


#{E^GBi) 


E 


def 

P 2 — 


1 


e G Er^nBi 


#(A„ n E 2 ) 


E 


eG E„nB2 


For all the examples we use the forcing terms h" = 0 for every n G N. The FEM approximation is done 
with 100 elements per edge with uniform grid. For each example we present two graphics for values of n 
chosen from {10, 20, 50,100, 500,1000}, based on optical neatness. For visual purposes in all the cases the 
edges are colored with red if K{e) = 1 or blue if K{e) = 2. Also, for displaying purposes, in the cases 
n G {10, 20} the edges ViVo are labeled with “f” for identification, however for n G {50,100, 500,1000} the 
labels were removed because they overload the image. 
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5.2. Verification Examples 

Example 1 (A Riemann Integrable Forcing Term). We begin our examples with the most familiar 

context as discussed in Remark 0 Define F : D — > IR by F(tcos£,tsui£) tt^ sin(7rt) cos(£). Since both 
sequences {vi : £ € N,£ = 0 mod 3} and {v{ : £ GTN,£^0 mod 3} are equidistributed, Weyl’s Theorem^ 
implies 

Fi = me[E|f2i,] = F 2 = me[E|Q^] = me[E] = 0. 


Here Fi, F 2 are the limits defined in Hypothesis For this case the exact solution of the upscaled 


Problem (26) is given by p p +p S Rq (^), with pi{t) = P 2 {t) = 0. For the diffusion coefficient we 
use the deterministic one, defined in (|33a[). The following table summarizes the convergence behavior. 


Example^: Convergence Table, K = Kd. 


n 

IIPi” -Pl||L2(eJ 

\\P2 -p2\\L^e,) 

ii^r - Pill 1 ( 62 ) 

IIP 2 " -P2 11 ^ 1 ( 62 ) 

10 

0.3526 

0.1717 

0.8232 

0.3216 

20 

0.0180 

0.0448 

0.0900 

0.0889 

100 

0.0160 

0.0059 

0.0395 

0.0116 

1000 

5.8352 X 10-4 

8.27772 X 10-4 

0.0012 

0.0016 



Figure 3: Solutions Example m Diffusion coefficient Kd-, see ( |33a[ >. The solutions depicted in figures (a) and (b) on the 
edges v^vq are colored with red ir Kd(viVQ) = 1 (i.e., £ = 0 mod 3), or blue if Kdiv^v^) = 2 (i.e., i ^ 0 mod 3). Forcing term 

F : Q, ^ IR, F{tcosi,tsini) sin(7rf) cos(.^). 


Example 2 (Probabilistic Flexibilities for Example^. This experiment follows the observations of 
i?emarfc[^ In this case take X Kp, defined in ( 33b| ). Let Z : N —>■ [—100,100] be a random variable with 

uniform distribution and define Y : £1g —> IR by F(tcos£, tsinf) tt^ sin(7rt) cos(£) + Z(£). It is direct to 
see that X and Y satisfy Hypothesis^ and due to the Law of Large Numbers, they also satisfy (31) and 
(321 respectively. Therefore 


Yi = me[E|f2i ] = Y2 = me[E|f22 ] = mg[F] = p^ = p^ = 0 . 
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The following table is the summary for a fixed realization of X (to keep the edge coloring consistent) and 
different realizations of Y on each stage. Convergence is observed, as expected it is slower than in the 
previous case. This would also occur for different realizations of X and Y simultaneously. 


Example^ Convergence Table, K = Kp. 


n 

\\Pl -plhHe^) 

\\P2 -p2\\L^(e£) 

\\Pi -PiWn^ie^) 

\\P2 -p2\\Hiie^) 

10 

0.5534 

0.0938 

1.6629 

0.5381 

20 

0.0965 

0.1594 

0.5186 

0.3761 

100 

0.0653 

0.1322 

0.3809 

0.2569 

1000 

0.0201 

0.0302 

0.0658 

0.0597 




Figure 4: Solutions of Example Fixed realization of diffusion coefficient Kp, see | |33b[ |. Forcing term Y : Qq —^ IR, 

Y(t cos t, t sin t) TT^ sin(7rt) cos(t) + Z{£), with Z : IN —> [—100,100], random variable and Z ~ uniformly. Different 

realizations for Y on each stage. The solutions depicted in figures (a) and (b) on the edges V£Vo are colored with red if 
Kp{vivo) = 1 (lEf/Cp = 1] = i), or blue colored if Kp(v£Vo) = 2 (lE[i^p = 2] = |). 


Example 3 (A non-Riemann Integrable Forcing Term). For our final theoretical example we use a 
non-Riemann Integrable forcing term. Moreover, the following function is highly oscillatory inside each 
subdomain and ^Iq, and it can not be seen as Riemann integrable when restricted to any of the sub- 
domains. Let F : ftc —f IR be defined by 


F{tcos£,tsmi) — 


sin(27rt) + (—1)L6J x 10 x (£ — ), 

TT^ sin(7rt) -h (-l)LsJ x 10 x (^ - [^J), 


£=0 mod 3, 
0 mod 3. 


(34) 


On one hand, both sequences {uf:£sN,£ = 0 mod 3} and {t'f:£sN,£^0 mod 3} are equidistributed. 
On the other hand both parts of the forcing term, the radial and the angular are Cesaro convergent on each 
£1q for i = 1,2. The Cesaro average of the angular summand is zero on for i = 1,2. In contrast, the 
radial summand can be seen as Riemann integrable separately on each for i = 1,2, therefore, due to 
Weyl’s r/ieorcTTi its Cesaro average is given by Fi = me[F|Q^] and F 2 = mg[F|Q^]; more explicitly, 

Fi(t) = (27r)^ sin(27rt) , P' 2 (t) = sin(7rt). 
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(35) 








For this case the exact solution p = plai +^10-2 G 1) of the upscaled Problem (26) is given by 


Pi{t) = sin(27rt), P2{t) = ^ sin(7rt). (36) 

We summarize the convergence behavior in the table below. 

Example^ Convergence Table, K = Kd- 


n 

Ibr -Pl||L2(ei) 

\\P2 -P2\\LHe,) 

ibr-Pi 1 + 1 ( 62 ) 

Ib2” -:P2|+i(e2) 

10 

1.6392 

0.4900 

5.7447 

1.3210 

20 

0.4127 

0.9305 

1.8930 

1.7782 

100 

0.2125 

0.3312 

0.4986 

0.6275 

1000 

0.0138 

0.0189 

0.0852 

0.0371 




Figure 5: Solutions of Example Diffusion coefficient K^, see l |33a[ l. The solutions depicted in figures (a) and (b) on the 
edges v^vo are colored with red if Kd{vivo) = 1 (i.e., £ = 0 mod 3), or blue if Kd{vevo) = 2 (i.e., £ ^ 0 mod 3). Forcing term 
F : flc —>■ IR see | |34^ . 


5.3. Numerical Experimentation Examples 

In this section we present two examples, breaking different hypotheses of those required in the theoretical 
analysis discussed above. As there is not a known exact solution, we follow Cauchy’s convergence criterion 
for the sequences {p" : n € IN} with t = 1,2. However, we do not sample only points but intervals of 
observation and report the averages of the observed data. More specihcally 


n 


1 

To 


n+5 

\\P^ - 


3=TL-i 


and 

for 




def 


1 

10 


n+5 

I] \\Pl 


j=n-4 


pr^UHe.), 


i = 1,2, n = 10,20,100,1000. 
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Example 4 (A Locally Unbounded Forcing Term). For our experiment we use a variation of Exampfe 
1^ keeping the well-behaved radial part but adding an unbounded angular part, which is known to be Cesaro 
convergent to zero. Consider the forcing term F : flo IR, defined by 


F(tcost', fsinf) 


47r^ sin(27rt)-I-(—1)^-\/£, £ = 0 mod 3, 

TT^ sm{Trt) + {—lyVi, £y0 mod 3. 


(37) 


Clearly, supgg^; ||E||^ 2 (g) = oo i.e.. Hypothesis is not satisfied. It is not hard to adjust the techniques 
presented in Section |4.1| to this case, when the forcing term is Cesaro convergent without satisfying the 
condition supgg^ ||F||i 2 (g) = oo; however the properties of edgewise uniform convergence of Section 
not be concluded. Consequently, we observe the following convergence behavior. 


3.2 


can 


Example^ Convergence Table, K = Kd- 


n 


n 

^2 

<5? 

^2 

10 

0.1547 

0.1578 

2.7336 

2.7338 

20 

0.0618 

0.0645 

1.0734 

1.0747 

100 

0.0277 

0.0224 

0.3394 

0.3320 

1000 

0.0086 

0.0065 

0.0984 

0.0955 



Figure 6: Solutions of Example [ 4 ] Diffusion coefficient K^, see l |33a| l. The solutions depicted in figures (a) and (b) on the 
edges V£Vo are colored with red if Kd{viVo) = 1 (i.e., £ = 0 mod 3), or blue if Kd{v£VQ) = 2 (i.e., £ ^ 0 mod 3). Forcing term 
F : Qg —> IR see | |37| l. 


Example 5 (A Forcing Term with Unbounded Frequency Modes). For our last experiment we use 
a variation of Example keeping it bounded, but introducing unbounded frequencies. Consider the forcing 
term F : D,c IR- defined by 


F{tcos£,tsm£) = 


47r^ sin(27rt • £), 
TT^ sin(7rt • £), 


£=0 mod 3, 
£y 0 mod 3. 


(38) 


Clearly, F verifies Hypothesis , then Lemma implies edgewise uniform convergence of the solutions, 
however Hypothesis-{u) is not satisfied. Therefore, we observe that the whole sequence is not Cauchy, 
although it has Cauchy subsequences as the following table shows. 
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Example^ Convergence Table, K = K^- 


n 

e? 

n 

^2 

<5? 

xn 

^2 

10 

0.0264 

0.0267 

0.4157 

0.3835 

20 

0.0078 

0.0089 

0.1342 

0.1327 

100 

0.0004 

0.0005 

0.0077 

0.0076 

500 

0.00004 

0.00004 

0.00073 

0.00072 

1000 

0.00066 

0.00049 

0.0081 

0.0078 

1200 

0.00004 

0.00005 

0.000787 

0.000786 




Figure 7: Solutions of Example [s] The Diffusion coefficient see ( |33a[ |. The solutions depicted in figures (a) and (b) on 
the edges v^vq are colored with red if K^iv^vo) = 1 (i.e., 1 = 0 mod 3), or blue if K^(v(Vo) = 2 (i.e., i ^ 0 mod 3). Forcing 
term F : Qq IR- see ||38|l. 


It follows that this system has more than one internal equilibrium. Consequently, an upscaled model of 
a system such as this, should contain uncertainty which, in this specific case, remains bounded due to the 
properties of the forcing term F. 

5.4- Closing Observations 

(i) The authors tried to find experimentally a rate of convergence using the well-know estimate 

log -pr\\- log \\pr-pr^w ■ 9 

* ^ log iipf -pr'll - log iiiPr' - 

The sampling was made on the intervals n — 5 < j < n -I- 5, for n = 10,20,100,500 and 1000. 
Experiments were run on all the examples except for Example In none of the cases, solid numerical 
evidence was detected that could suggest an order of convergence for the phenomenon. 

(ii) Experiments for random variations of the examples above were also executed, under the hypothesis 
that random variables were subject to the Law of Large Numbers. Convergence, slower than its 
corresponding deterministic version was observed, as expected. This is important for its applicability 
to upscaling networks derived from game theory, see |20) . 
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6. Conclusions and Final Discussion 


The present work yields several accomplishments and also limitations as we point out below. 


(i) The method presented in this paper can be easily extended to general scale-free networks in a very 
simple way. First identify the communication kernel (see |13]b Second, for each node in the kernel, 
replace its numerous incident low-degree nodes by the upscaled nodes together with the homogenized 
diffusion coefficients and forcing terms, see Figure 

(ii) The particular scale-free network treated in the paper i.e., the star metric graph, arises naturally in 
some important examples. These come from the theory of the strategic network formation, where the 
agents choose their connections following utilitarian behavior. Under certain conditions for the benefit- 
cost relation affecting the actors when establishing links with other agents, the asymptotic network is 
star-shaped (see ED). 

(iii) The scale-free networks are frequent in many real world examples as already mentioned. It follows that 
the method is applicable to a wide range of cases. However, important networks can not be treated the 
same way for homogenization, even if they share some important properties of communication. The 
small-world networks constitute an example since they are highly clustered, this feature contradicts 
the power-law degree distribution hypothesis. See [55] for a detailed exposition on the matter. 

(iv) The upscaling of the diffusion phenomenon is done in a hybrid fashion. On one hand, the diffusion on 
the low-degree nodes is modeled by the weak variational form of the differential operators defined over 
the graph, but ignoring its combinatorial structure. On the other hand, the diffusion on the communi¬ 
cation kernel will still depend on both, the differential operators and the combinatorial structure. This 
is an important achievement, because it is consistent with the nature of available data for the analysis 
of real world networks. Typically, the data for central (or highly connected) agents are more reliable 
than data for marginal (or low degree) agents. 

(v) The central Cesaro convergence hypotheses for data behavior (stated in Lemma(iii), as well as those 
contained in Hypothesis in order to conclude convergence have probabilistic-statistical nature. This 
is one of the main accomplishments of the work, because the hypotheses are mild and adjust to realistic 
scenarios; unlike strong hypotheses of topological nature such as periodicity, continuity, differentiability 
or even Riemann-integrability of the forcing terms (see m)- This fact is further illustrated in Example 
[^ where good asymptotic behavior is observed for a forcing term which is nowhere continuous on the 
domain Hg of analysis. 


(vi) An important and desirable consequence of the data hypotheses adopted, is that the method can 
be extended to more general scenarios, as mentioned in Remark reported in Subsection 5.4 and 
illustrated in Examples Hi and i Moreover, Example i suggests a probabilistic upscaled model for 
the communication kernel, to be explored in future work. 


(vii) A different line of future research consists in the analysis of the same phenomenon, but using the 
mixed-mixed variational formulation introduced in |23j instead of the direct one used in the present 
analysis. The key motivation in doing so, is that the mixed-mixed formulation is capable of modeling 
more general exchange conditions than those handled by the direct variational formulation and by the 
classic mixed formulations. This advantage can broaden in a significant way the spectrum of real-world 
networks which can be successfully modeled and upscaled. 


(viii) Finally, the preexistent literature typically analyses the asymptotic behavior of diffusion in complex 
networks, starting from fully discrete models (e.g., |M1[1S|). The pseudo-discrete treatment that we 
have followed here, constitutes more a complementary than an alternative approach. Depending on 
the availability of data and/or sampling, as well as the scale of interest for a particular problem, it is 
natural to consider a “blending” of both techniques. 
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